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TECHNICAL NOTE 3341 


AN ANALYTICAL ESTIMATION OF THE EFFECT OF TRANSPIRATION 
COOLING ON THE HEAT-TRANSFER AND SKIN-FRICTION 
CHARACTERISTICS OF А COMPRESSIBLE, 

TURBULENT BOUNDARY LAYER 


By Morris W. Rubesin 
SUMMARY 


An analysis based on mixing-length theory is presented which 
indicates that surface blowing associated with transpiration cooling 
systems produces large reductions in both the heat-transfer and skin- 
friction coefficients Рог a turbulent boundary layer on a flat plate. 
The numerical results are restricted to the case of air blowing into 
air. Тһе effects of blowing are indicated to be similar for high-speed, 
compressible flow to those for low-speed, incompressible flow. 


INTRODUCTION 


The frictional heating of the outer surfaces of high-speed aircraft 
has become a major problem in the design of these aircraft. Without 
thermal protection and at equilibrium conditions, the surfaces of these 
aircraft will begin reaching intolerable temperatures at Mach numbers 
even as low as 2. The designer, therefore, must provide thermal 
protection for his aircraft's skin. This protection can be achieved in 
several ways; for example, by altering the aircraft's shape to avoid 
sharp or pointed frontal surfaces, by providing a cooling system for the 
skin, and by providing a protective thermal insulating layer between 
the hot air in the boundary layer and the skin. A transpiration cooling 
system, one in which the coolant passes through small pores in the skin 
and into the outside boundary layer, has the advantage of providing both 
cooling of the skin and a protective insulating layer of coolant. It 
appears that these advantages make a transpiration cooling system most 
effective (ref. 1}. It is noted that the advantage of evaporation can 
also be incorporated into a transpiration cooling system. 


Much of the available work dealing with transpiration cooling is 
restricted to analyses dealing with the laminar boundary layer. The 
literature is quite extensive and references 2 and 3 represent examples 
of these analyses. Investigations which are more allied to the work 
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presented in this paper are represented by references 4, 5, and б.” 
These analyses treat the case of the turbulent boundary layer by divid- 
ing the boundary layer into two parts, а laminar sublayer and а fully 
turbulent outer region. The end results of these analyses, in effect, 
relate the coefficient of heat transfer under conditions of transpira- 
tion with two parameters; namely, (1) the velocity &t the interface of 
the sublayer and outer turbulent portion and (2) a Reynolds number based 
on the distance of this interface from the surface with properties 
evaluated either at the wall or free-stream temperature. The rate of 
transpiration acts as an independent variable in the latter relation. 
Because it is known that both of these parameters are dependent on the 
local skin friction for the case of zero transpiration, it would be 
expected that their partial dependence on the local skin friction would 
continue even under the conditions of transpiration. Іп effect, then, 
the results of these analyses are modified Reynolds analogies relating 
heat transfer and skin friction when blowing occurs. A similar analysis 
comprises a portion of this report; however, it is modified to also 
include the frictional dissipation occurring in the high-speed air flow 
over the surfaces of aircraft. 


Several experiments have been performed concerning transpiration 
cooling in tubes or channels having turbulent boundary layers (refs. 5, 
7, and 8). Although the results point out the advantages of transpira- 
tion cooling, the geometry of the tests and the low velocities employed 
in the boundary layer make the results of these tests too specific to 
apply in the general case of aircraft. The experiment described in 
reference 9, although also employing low air speeds, provides data on 
а flat plate, which is а fundamental aerodynamic shape. These data will 
be compared with the present analysis in & later portion of this report. 


From this brief review of literature dealing with transpiration 
cooling, it is apparent that there is no information about the influence 
of transpiration on the behavior of а compressible, turbulent boundary 
layer, such as exists on the surface of a high-speed aircraft. The 
purpose of the present report is to present an approximate analysis for 
determining the effect of transpiration on a high-speed turbulent 
boundary layer. Because of the many uncertainties inherent in the 
analysis, it mainly has heuristic value. Experimental data will be 
required before the value of the analysis as a means of interpolating 
or extrapolating limited amounts of data can be assessed. 


+A simultaneous investigation, essentially the same as is contained 
in portions of the present report, has been reported recently in an 
article by W. H. Dorrance and F. J. Dore, entitled "The Effect of Mass 
Transfer on the Compressible Turbulent Boundary Layer Skin Friction and 
Heat Transfer." Jour. Aero. Sci., vol. 21, no. 6, June 195}. 
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SYMBOLS 


arbitrary constant 


parameter, 


arbitrary constant 
1 + 7-3 м. 2 
parameter, ес 
Ty/Tı 


specific heat at constant pressure 


local skin-friction coefficient 


: St 2 Е 
ratio defined by equation (10 
3 сг/Р 5? у €q ( ) 


elliptic integral of the first kind 


exponential term defined by equation (A33) and (АЗ!) in 
Appendix А 


function defined by equation (6) 
local heat-transfer coefficient defined by equation (А78) 
thermal conductivity 


parameters of elliptic integral of first kind (see eqs. (A66) 
and (А68)) 


mixing length parameter, 1 = Ky 
mixing length 


Mach number 


e 
Prandtl number, = 


| € 
turbulent Prandtl number, A 


Rg 


St 


Pins Pn 
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local heat-transfer rate per unit area 
local temperature recovery factor (see eq. (A80)) 


Reynolds number based on distance along plate and free-stream 
fluid properties 


Reynolds number based on momentum thickness and free-stream 
fluid properties 


Local Stanton number, 


Озу есуу 
temperature 
velocity component parallel to plate 
velocity component perpendicular to surface of plate 
coordinate parallel to axis of plate 
coordinate perpendicular to surface of plate 
parameter defined by equation (A58) 
parameter defined by equation (A59) 
ratio of specific heats, 1.40 for air 
eddy thermal conductivity, defined by equation (A11) 
eddy viscosity, defined by equation (A10) 


constant relating transpiration rate with local skin-friction 


coefficient, Py Vy 01x 


сг/2 
momentum thickness ш 
viscosity 
mass density 
parameter defined by equation (A60), 67> 
local shear stress 


parameter of elliptic integral of first kind 
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а) exponent in viscosity - temperature relationship 
Superscripts 
instantaneous fluctuating quantity 
- quantity &veraged with respect to time 
4 dimensionless quantity (see eqs. (А26) and (А27)) 
м dimensionless quantity (see eq. (A19)) 
Subscripts 
a interface of laminar sublayer and outer turbulent region 
c coolant fluid initial condition 
142 
е temperature recovery condition, Те = T, + r =>- 
2 Cp 
W surface 
1 outer edge of boundary layer 
ANALYSIS 


i Restrictions on Analysis 


The present analysis has several restrictions which have been 
imposed because of inherent difficulties or for purposes of simplifica- 
tion. These restrictions are as follows: 


l. Only the compressible, turbulent boundary layer is considered 
(usual boundary-layer assumptions). 


2. The boundary layer is on a flat plate placed parallel to the 
free=stream direction. 


3. The coolant fluid is the same as the boundary-layer fluid and 
enters the boundary layer at the temperature of the surface 
of the plate (no diffusion). 
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4. Тһе boundary layer is idealized to be composed of two distinct 
sections: a laminar sublayer next to the surface, and а 
fully turbulent outer layer. (Note: It is assumed that 
velocity, temperature, heat flux, and shear are continuous 
at the interface of the two layers.) 


2. In the fully turbulent outer layer, the transport of heat is 
proportional to the transport of momentum. 


Parts of the Analysis 


For relative simplicity the analysis is divided into two parts. 
The first part contains a determination through the use of mixing-length 
theory of the effect of transpiration on skin friction and heat transfer 
іп а compressible turbulent boundary layer when the Prandtl number is 
unity. The assumption Pr = 1 greatly simplifies this portion of the 
analysis. The second part contains a determination of the effect of 
transpiration on the relation between the coefficients of skin friction 
and heat transfer when the Prandtl number is not unity. An expression 
for temperature recovery factor is also found in the second part. 


A combination of the results of the two parts (the product of ce 
when Pr = 1 and St/cp when Pr £ 1) allows the determination of the local 
heat-transfer coefficients under the conditions of transpiration when 
Prandtl number is not unity. То permit this combination, it is necessary 
to employ the premise that the local skin-friction coefficient is 
relatively independent of Prandtl number. For the case with no tran- 
spiration, this premise is good Рог a Prandtl number which corresponds 
to air (ref. 10). It is not expected that small amounts of transpiration 
should &lter the basic premise. 


Description of Method of Analysis 


This description of the analysis is intended to point out its 
salient features without introducing the complexities of the detailed 
mathematics. The mathematical details may be found in Appendix A. 


The basic equations that are used in this analysis are the Navier- 
Stokes equations representing conservation of momentum, the conservation 
of energy equation, and the continuity or conservation of mass equation. 
When all the quantities in these equations are represented as the sum 
of a time averaged quantity plus an instantaneous quantity, these 
equations can be made to apply to turbulent flow (ref. 11). When these 
"turbulent flow" equations are time averaged and the usual boundary- 
layer assumptions are imposed, the resulting boundary-layer equations 
are (ref. 10): 
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Barred quantities represent time averages, while primed quantities 
represent instantaneous values of fluctuating quantities. While the 
specific heat is maintained constant, the fluid properties of density, 
viscosity, and thermal conductivity are allowed to vary with temperature. 
Frictional dissipation is also included. 


Equations (1) through (3) are intractable at present, because of 
their complexity and because the behavior of the turbulent fluctuating 
components in shear flow is not known. Certain nonrigorous simplifica- 
tions are necessary to achieve a solution useful from an engineering 
viewpoint. | 


The first simplification necessary in achieving an engineering 
solution to the problem is to assume in equations (1) through (3) that 
the variations of the dependent variables with respect to x are 
negligible compared to their variations with respect to y in the 
evaluation of the local velocity distribution. This allows replacing 
partial differential equations by much more easily solved ordinary 
differential equations, with y ав the independent variable. The above 
assumption is based on experience with the low-speed, incompressible 
skin-friction problem without transpiration. It has been found that 
this assumption combined with a group of assumptions concerning the 
mechanism of turbulence and the structure of the boundary layer together 
with certain mathematical simplifications can be made to yield skin- 
friction coefficients which agree well with the experimental results. 


At this point it is necessary to define the Prandtl number as it 
is used in this report. Іп the conventional definition Prandtl number 
is the ratio of kinematic viscosity to the thermal diffusivity of a fluid 
where the transport of momentum and heat is promoted by molecular means. 
This Prandtl number is then a controlling variable for processes con- 
taining laminar flow where molecular transport phenomena occur. When 
the flow is turbulent, however, the large-scale transport of momentum 
and heat is promoted by the eddying motion, and the conventional Prandtl 
number does not govern this mechanism. In order to relate the turbulent 
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transport of momentum and heat transfer, it is necessary to define а new 
term which has the same characteristics as Prandtl number, except that 
properties defining the term are based on eddying motion. To distinguish 
this term from the conventional Prandtl number, it is called the turbulent 
Prandtl number. А turbulent Prandtl number of unity is equivalent to 

the Reynolds analogy. 


When the basic equations are reduced to ordinary differential 
equations, they are applied to the fully turbulent outer layer by 
omitting all the terms which are based on molecular transport. In the 
laminar sublayer, however, all the terms based on eddying transport are 
omitted. Thus, two sets of equations are employed, the solutions of 
each being matched at the interface of the two layers. Experience has 
shown that this two-layer model of a turbulent boundary layer is quite 
good as long as both the Prandtl number and turbulent Prandtl number are 
close to unity. 


At this stage, the analysis is divided into two parts. In the 
first part it is assumed that Prandtl number and turbulent Prandtl number 
are both equal to unity, which results in a direct relationship between 
the temperature and velocity. Because of this relationship, it is not 
necessary to solve the energy equation. The momentum equation is solved 
together with the continuity equation to yield the velocity (and tem- 
perature) distribution through the boundary layer with the local surface 
shear as & parameter. The velocity distribution is then substituted in 
the integral expression for the momentum thickness. As this integral 
is not readily solvable in closed form, the first term of a series 
solution is employed. The first term represents an asymptotic solution 
for the case of extremely low values of the skin-friction coefficient 


and low rates of transpiration; that is, / cp/2 << 1 and 


J (ey vy/0131) << l. Thus, & relationship is obtained between the local 


skin-friction coefficient and Reynolds number with the momentum thickness 
ав the characteristic dimension. When the von Karman momentum-integral 
equation is integrated using the relationship се(Вр), there results an 


expression for the local skin-friction coefficient in terms of а Reynolds 
number using the length along the plate as the characteristic dimension. 
This latter integration is also not made in closed form, the first term 
of a series expansion again being used as an asymptotic relationship 

for small values of the skin-friction coefficient and the transpiration 
rate. The first part of the analysis, consequently, results in the 
determination of the effect of transpiration on the relationship 

сг(Вх} when Pr and Pr, are both unity. 


It is known that the Local heat-transfer coefficient is dependent 
on Pr and Pry. To evaluate this dependence, the second part of the 
analysis is performed to determine the relationship between the local 
Stanton number and the local skin-friction coefficient under conditions 
of transpiration when Pr and Pr; are not unity. This is achieved by 
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taking the ratio of members of the energy equation to corresponding 
members of the momentum equation. The ratio of the members forms an 
ordinary differential equation relating total energy and velocity. When 
the differential equation is integrated in closed form across the bound- 
ary layer, there is obtained a relationship between Stanton number and 
the local skin-friction coefficient. An ancillary result is an expres- 
sion for the local temperature recovery factor. 


A combination of the two parts of the analysis provides St(R,) 
under conditions of transpiration for Prandtl numbers other than unity. 


Summary of Results of Analysis 


The results of part I of the analysis can be summarized as 
follows: 


RoK 
0 (4) 
PrrVtr/P1U 
yt 1+ W 37 1-1 
8. ср/2 


ОГ 


(5) 


where 
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А ам: В = Сү, 
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when 


or (6) 


when 


Note: 
B- VB? + ААХ 


oA? 
„вам B^ + ЦАР 
QA? 
7 (сг/2) 


РУ 
ра Ча. 


The term F(k,?) is an elliptic integral of the first kind. 


For the values of К, уі, and Ug to be used in the above see "DISCUSSION" (note Ùa = па 2/2). 


The results of part II of the analysis can be summarized as follows: 


РУ и: 
білі 2 


wy ) 
St Ж PLUL | ( T) 
Cp l1-Pr Pr-Pr. 
Mun ( Ср DwYw "T Ср ) Ср 
2 PLU 2 2 


(Pr,-Pr) 
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(8) 
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DISCUSSION 


Evaluation of К, ug, Ye, 


The end results of the analysis, as represented by equations (4) 
through (8), are entirely general until specific values of K, ug, and 
ye аге selected. This means that before numerical answers can be 
determined from this analysis, information other than that contained in 
the analysis must be employed. This information, that is, the values 
of К, ug, and уд, must be determined from experiment. Thus, the 
numerical results of this report can at best be as good ав the experi- 
ments used in selecting К, ui, and У. 


Low-speed case with no transpiration.- For the incompressible, zero 
transpiration case К, ug, and ух can be evaluated in three ways. The 
first way is from velocity-distribution data, where К is а measure of 
the rate of increase of the velocity with distance from the surface, 
and ug and уҙ are the conditions at the interface of the sublayer and 
outer turbulent portion. The second way is from cp(Rg) data. In this 
evaluation the terms К, 45, ала Yà lose some of their physical signi- 
ficance in that they may be considered as arbitrary constants which 
absorb some of the inaccuracies introduced by the mathematical simpli- 
fications employed in evaluating the momentum thickness from the 
velocity-distribution data. The third way is from ce(Ry) data where, 
again, К, ug, and уі are treated as arbitrary constants, absorbing 
still other mathematical simplifications. 


The experimental velocity-distribution data can be plotted in the 
dimensionless form ut vs. у’ when local surface shear stresses have 
been obtained simultaneously. By choosing К, ug, and yl appropriately, 
the analysis can be made to fit portions of the boundary-layer data very 
well. If the quantities K, Чи» and Ya are adjusted so that the inner 
portion of the turbulent layer is well represented, the analysis, how- 
ever, will not fit the data in the outer portion of the layer, and vice 
versa. This lack of agreement between analysis and experiment is not 
entirely surprising when it is recalled that the analysis is based on 
& primitive mixing-length theory and that several mathematical simpli- 
fications have been imposed.  ' 


Before the experimental data of ср(Во) ог ер(Вх) сап be used to 
evaluate the arbitrary quantities, it is necessary to relate ug and yi 
analytically. This is necessary because the skin-friction data will 
only supply either K and Us or K and Ул with sufficient accuracy. 
For the incompressible, no transpiration case, & comparison of the 
quantities evaluated in the three different manners is shown in the 
following table: 
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From velocity 
survey data 

(inner portion of 

turbulent layer) 


From cp(Rg) | From сг(В,) 


These values were evaluated from figures of data summarized in reference 
12. It is apparent that the values of K and u$ given by the alterna- 
tive methods do not agree, the extreme values of К differing by about 
14 percent and the extreme values of uà also differing by about 14 per- 
cent. These differences, undoubtedly, are due to the inadequacies of tbe 
mixing-length theory and the mathematical approximations imposed. 


For the purposes of the present report, however, tbe chief concern 
is adequate representation by the analysis of cg(Ry4). In figure 1, 
there has been reproduced & figure from reference 12 showing available 
data of с (Ех). The curves given by the present analysis are represented 
by solid lines. The analytical line based on cg(Rg) data differs from 


the c (Ry) data by less than 7 percent. The analytical line, based оп 


velocity-distribution data, is about 12-percent higher than the c (Rx) 
data. From this it can be inferred that no matter which method is 
employed to obtain К, ug, and yi, use of these values to obtain ec p(R,) 


will result in answers adequate for engineering purposes. 


Low-speed case with transpir&tion.- There exist insufficient data 
to define precisely the effect of transpiration on К, ui, and уг The 


data of reference 9 do not allow these determinations even though they 
include temperature=distribution surveys, velocity-distribution surveys, 
and measurements of local heat transfer. The temperature-distribution 
data were not used because evaluating K, оў, апа Ya from temper&ture- 
distribution surveys and local surface heat-transfer measurements 
requires making assumptions concerning the analogy between heat transfer 
and momentum transfer. Examination of the velocity-distribution data 
revealed that local skin-friction coefficients were more than 20 percent 
higher than the values in figure 2 for the case with zero transpiration. 
These high values were attributed to surface roughness effects. Since 
it is not known how roughness is affected by transpiration, it is 
believed that the velocity-distribution data may include the effects of 
both. Because heat transfer in turbulent flow seems to be affected by 
roughness to a much lesser extent than is skin friction (ref. 13), it 
was decided to employ the heat-transfer data in the evaluation of K, 
us and Vers As the data are not sufficiently precise to do this 
directly, it was necessary to employ the alternative approach of postu- 
lating values of K, Uns and У and then of comparing the end results 
with the data of reference 9. Physical considerations were employed in 
guiding these postulates. 
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From the physical viewpoint it would be expected that Ug and. Ya 


should be consistent with the velocity-distribution curve within the 
sublayer, equation (A25). In addition it would be expected that the 
Reynolds number at the outer edge of the sublayer, based on the thick- 
ness of the sublayer and equal to оауё, should be nonincreasing under 
the disturbing influence of transpiration. Consequently, the two 
assumptions investigated are: | 


(&) Uayt = (13.1)? solved simultaneously with equation (A25) and 
(b) ug = 13.1 and E from equation (A25). 


The first of these represents the case of а constant sublayer Reynolds 
number, while the second corresponds to & Reynolds number which becomes 
smaller as the transpiration rate is increased. A comparison of the 
individual results based on these &ssumptions provides ап indication of 
the influence the assumptions concerning tbe sublayer may have. То 
examine this point further, other numerical results were obtained using 
the condition that ug = Yà = 13.1, identical with the no transpiration 
cage. It should be noted that this latter condition is inconsistent 
with the sublayer equation (eq. (A25)). 


Besides the question of the boundary conditions at the outer edge 
of the sublayer, there exists the problem of choosing a value of K, s 
the mixing-length parameter, under the conditions of transpiration. Only 
well-defined experiments can &nswer this problem and, &t present, the 
only recourse available is to assume it is unchanged by the transpiration; 
namely, К = 0.392. In addition to these assumptions, the conditions 
M = 0, T/T, = l, Pr = 0.72, and Pry = 1.00 are imposed on the analysis 
to make it conform with the test conditions of reference 9. 


Comparison with Data’ 


A comparison of the data and analysis is shown in figure 2 for the . 
case where the transpiration rate, бүтүр / Оза. › is equal to zero. The 


local Stanton number is plotted as a function of Reynolds number with 
the distance along the channel wall used &s the characteristic dimension. 
It appears that the data of run H-1 correspond to & laminar boundary 
layer; the data of runs H-2 and H-3 correspond to а boundary layer which 
is laminar at the lower Reynolds numbers and transforms to one that is 
turbulent at the higher Reynolds numbers; and the data of run Н-! 
correspond to a boundary layer which has become turbulent very close 

to the leading edge of the channel wall. When the data of run Hel are 
compared with Pohlhausen's theoretical result for laminar flow, it is 
noted that the data are generally about 20-percent higher than theory. 
Since Pohlhausen's theory was checked experimentally in the past by Fage 
and Falkner (ref. 14}, some doubt is raised concerning the accuracy of 
the data with laminar flow. The data corresponding to turbulent flow, 
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however, agree fairly well with the well-known empirical Colburn curve, 
reference 15, or with a curve derived from the present analysis (note 
for case where py Vy/ 013 = 0, various postulates about the sublayer are 


identical). Тһе deviations shown by the high Reynolds number data for 
each of runs H-3 and Н-Ё аге not considered. 


For the purpose of the present paper, the main point is considera- 
tion of the agreement of the data corresponding to turbulent flow with 
the results of the present analysis. It should be emphasized that for 
comparison purposes the data should actually be plotted against a Reynolds. 
number in which the characteristic dimension is the effective length. of 
& fully turbulent boundary layer. The effective length is defined as 
the distance &long the flat plate where а boundary layer, fully turbulent 
from the start, would have the same Stanton number as the actual boundary 
layer for which transition from laminar to turbulent flow has occurred 
downstream of the leading edge. The data given in reference 9 sre 
insufficient to define precisely the effective start of turbulence, 80 
it is necessary to resort to estimates. Because the Stanton number 
approaches infinity as Reynolds number approaches zero, any of the data 
points corresponding to turbulent flow having a finite value of Stanton 
number are downstream of the effective starting point of the turbulent 
boundary layer. Thus, if the effective starting point of a fully tur- 
bulent boundary layer is chosen as the point where the data indicate 
transition has just ended and the flow is completely turbulent, the 
effective-length Reynolds number so determined would certainly be smaller 
than the true effective Reynolds number. When this extreme correction 
is applied to the data, that is, 1.8х107 is subtracted from the Ry | 
for run H-2, 9Х105 is subtracted from the Ry for run H-3, and the Ry 
of run H-4+ are left unchanged, it is found that the data cluster near 
the analytical result with a spread of about 20 percent as shown in 
figure 3. The results based on the correct effective starting length 
would undoubtedly lie between those shown in figures 2 and 3 and would 
be in good agreement with the proposed analytical method, considering 
the inherent scatter in the data. 


Before & comparison of the proposed analytical method and the data 
can be made under conditions of transpiration, it is necessary to discuss 
aspects of the variation of the transpiration rate along the surface of a 
plate. If it is assumed that both conduction within the porous surface 
and radiation are negligible, a heat balance on & portion of the porous 
surface results in the following equation: 


PwVW 


piu; 


(Ту - Те) (9) 


St(Te - M) = 


Because St varies with х, it is apparent that if Te and Te are 
constant with respect to distance along the plate, the only condition 
compatible with & constant surface temperature is DwVw/DPiUui = 65% 
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where ес is а constant of proportionality. On the contrary, if 
pu vyu/o 101 = constant, the surface temperature Ty will vary along the 


surface of the plate. 


In the present analysis it is quite inconvenient to treat the first 
of these cases exactly. This arises from the fact that the analysis is 
divided into two parts, that is, the local Stanton number is determined 
from a local skin-friction coefficient that was determined for the 
particular case of Pr = Pr, = 1 and а prescribed variation of 
ouvy/o 111. Thus, to make бү /О iux = eSt would require ап iteration 


process. For convenience, and because it does not produce а large 
error in the practical case, the relationship 


Pr, Pra, ОҒЫ м) (10) 
Quem 


PLUL 
will be assumed. This assumption implies that the Reynolds number 
dependence of the function f is negligible. It will be shown later 
that although this is not exactly true, the assumption only introduces 
an error of about tl percent over the Reynolds number range of the 
experiments in reference 9. When equation (10) is introduced into 


equation (9) there results: 
ШЕ: | 
(Ty - Те) (11) 
CD. 
Thus letting Ои /о 11 = 6(сғ/2) results in essentially a constant 
surface temperature, where б is а constant of proportionality. 


For the case where the fluid properties are essentially constant, 
that is, low-speed flow with small temperature variations, it is possible 
to apply the analysis directly to the case where P гн /о 101 = constant. 


Some error will occur because of neglecting the effect of ах1а1 
variations in the surface temperature on the convective heat transfer. 
For the case with zero transpiration, it is known that continuous 
variations of surface temperature have little effect on heat transfer 
in а turbulent boundary layer (ref. 16). It is not expected that small 
amounts of transpiration should markedly alter this. Thus, the results 
provided by the analysis for Өрү / Рала = constant should be applicable 


for comparison with the data of reference 9. 
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Comparisons of the analytical results with the data of reference 9 
are presented in figures 4, 5, б, and T. Figures № and 5 correspond to 
constant oy vy/0 ааа, and figures 6 and T correspond to 


Ou y/o 102 = €(cp/2). In these figures no corrections for the effective 
starting length have been made for the Reynolds numbers of the data. 


In figure 4 the variation of Stanton number with length Reynolds 
number is shown for Py Vw/01ui = 0.002. The analytical results based 


on the assumed conditions of К, ug, and yf are also included. It is 
surprising how well either postulate compatible with the sublayer 
equation compares with the data. The analysis with either of these 
postulates gives about the correct order of magnitude for the reduction 
in Stanton number due to transpiration and, in addition, indicates the 
larger variation in Stanton number with changes in Reynolds number. The 
differences in the analytical results and the data are shown much more 
clearly in figure 5 where Ө үрт / Ола. = 0.006. In this figure neither 


analytical result based on either postulate agrees well with the data, 
however, both analytical results indicate a large reduction in St 

and a larger variation of St with Ну, as is exhibited by the data. It 
is interesting that the two analytical results bracket the data. In both 
figures, however, the postulate where the sublayer is considered unchanged 
by transpiration leads to analytical results which exhibit the poorest 
agreement with the data. Another point of interest is a confirmation 
that variations in surface temperatures affect the data to a small 
degree. This is seen from а comparison in figure № of runs 17b and 20b 
with runs 16а, 178, апа 20а, where the differences between а constant 
surface temperature and а variable surface temperature are masked by 

the scatter of the data. The results of figure 5 further substantiate 
this point. 


In figures 6 and 7 the variation of Stanton number with length 
Reynolds number is shown for РУМО 13 = 6(сғ/2). In view of the 


previous results, only the analytical results based on postulates 
compatible with the sublayer are compared with the data. To agree with 
the experimental values of € «їе values of. € for the two cases are 
1.04 and 1.15, respectively, in figure 6. Тһе corresponding values of 
Е in figure T are б = 2.68 and 3.34. Although the data scatter 
markedly, a comparison of the data and analysis is quite favorable when 
consideration is made of the probable effective starting point of the 
data of run Н-10. 


From the above comparison of data of reference 9 with the results 
of the present analysis, it can be concluded that for the low speed, 
incompressible case, answers sufficiently accurate for engineering 
purposes can be provided by the analysis when the arbitrary constants 
of the analysis are set equal to: 
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0.392 
ца = 13.1 


and y, is determined from equation (A25). It will be noted that a 
postulate considering the sublayer Reynolds number variable is employed 
in preference to a postulate where it is considered constant. This is 
done from the viewpoint of providing a more conservative answer (i.e., 
less reduction in St with transpiration) than is exhibited by the bulk 
‘of the data. 


Examples of Results of Analysis, Including Extensions 
to High-Speed, Compressible Flow 


Examples of the effect of transpiration on the Local skin-friction 
coefficient, the ratio of Stanton number to the skin-friction coefficient, 
and the local recovery factor are shown in figures 8, 9, and 10. These 
curves include cases of Mach number equal to О, 2, and i. Тһе curves 
shown are subject to the following conditions: 


(а) К = 0.392 


T 
(b) ut = 13.1 = 


1 
(c) yg is determined from equation (425) 
(4) ш=о.8 

For figure 8 the conditions for Prandtl number are: 
(e) Pr = 1.00 
(f) Pr, = 1.00 

In figures 9 and 10 the conditions for Prandtl number are: 
(в) Pr = 0.172. 
(а) Pr, = 1.00 


In all the figures, the parameter of the curves is the dimensionless 
B (owvw/01va) 


(ер/2) 


implying a varying transpiration rate along the surface, was chosen for 
the examples from the practical viewpoint that it requires a smaller 
amount of transpiration air to maintain a certain maximum surface 


transpiration parameter 6 . This condition of constant 6, 
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temperature over а region than does a constant transpiration rate 
system. Two surface-temperature conditions are presented: the first, 
is where the surface is at the total temperature and the second is where 
the surface temperature is at the free-stream static temperature. These 
surface temperatures are certainly the extreme limits of the range of 
temperatures likely to be encountered on the surface of a cooled aircraft. 


Conditions (а) through (с) are essentially the same as postulate (b) 
(page 14) defined previously, except that ug is based on fluid prop- 
erties evaluated at the surface temperature. This accounts for the 
temperature-ratio term in (b). Condition (d) represents the viscosity 
exponent corresponding approximately to conditions occurring in flight. 
Conditions (e) and (f) were chosen for purposes of mathematical simpli- 
fication while (g) corresponds to air and (h) implies an exact Reynolds 
analogy in the outer turbulent layer. 


It is apparent from figure 8 that transpiration has а very strong 
effect on lowering the value of the local skin-friction coefficient at 
&ll the Mach numbers considered &nd эф both surf&ce-temperature con- 
ditions. Ав an example, at M = 4 where Ту = Ті, the local skin-friction 


coefficient at Вх = 107 is reduced to around 15 percent of its original 
value by a blowing rate of pyvy/o1ur = 0.003 (Е = 20). 


From figure 9 it is seen that the effect of transpiration at all 
the conditions considered is to increase the ratio of Stanton number 
to the local skin-friction coefficient when Pr = 1. This increase is 
of a much smaller magnitude than the decrease in the local skin-friction 
coefficient, with the net result being that the Stanton number is also 
markedly reduced by transpiration. For example, at the point previously 
considered, the Stanton number falls to about 21 percent of its original 
value. 


In figure 10 there is shown & series of curves for the local 
recovery factor аб М = О. This figure is representative of all the 
other conditions considered, and the following discussion applies in 
general. From the figure it is noted that transpiration also affects 
the recovery factor; however, it can be observed that the effect of 
transpiration is not too large. Іп fact, the transpiration effect is 
much less than'the apparent effect of Reynolds number for the zero 
transpiration case. Since it is known from experiment (ref. 10) that 
the recovery factor at zero transpiration is relatively independent of 
Reynolds number, doubt is shed on the validity of the recovery-factor 
determinations of this analysis. From a consideration of relative 
results, all that can be concluded is that transpiration probably does 
not affect the recovery factor by a large amount and that a value of 
г = 0.9 should suffice in design calculations. 
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In applying the results of this analysis as expressed in figures 8 
and 9, it is convenient to express the parameter б in terms of the 
temperatures controlling a given cooling problem. If an effectiveness 
of cooling is defined by 


Ta - T 
ко 72. W (12) 
Te - Te 


for a surface where a heat balance (eq. (9)) applies, it is possible to 
relate б апа E by 


Е = — (13) 
С + 
(cp/2) 
Thus, once E is specified for a given design, it is possible to 
determine € and, consequently, cg, St, and руту for given flight 


conditions and surface temperature. 
CONCLUDING REMARKS 


It сап be concluded from the results of the present analysis that 
a transpiration cooling system is highly effective, even at high flight 
speeds. Not only does the coolant absorb & maximum amount of heat by 
attaining as its end temperature the temperature of the surface that it 
is cooling, but it has the further advantage that the effect of the 
transpiration greatly reduces both the skin friction and the amount of 
heat entering the surface. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 12, 1954 
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APPENDIX A 
DETAILS OF ANALYSIS 


When it is assumed that changes in the x direction are negligible 
compared to changes in the y direction, equations (1) through (3) 
become 


a (09 + рап) = 0 (A1) 
(6v + pv’) Mit yo m (A2) 
dy dy 
(бу + p'v!) i (ат. )-2 ау ES + 3 боро + араз) 
(A3) 


These equations can be simplified by integrating equation (Al). Тһе 
integr&l of equation (A1) is 


pv + p'v! = constant (Ab) 


At the surface of the plate 


ру = Dg Vw (A5) 
and 
p'v! = 0 (А6) 
Therefore, 
DV + piv’ = руту (АТ) 


Equations (A2) and (АЗ) thus become 


21 - би A8 
они ay ay AC ay © e 
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а в Ле. а Исав та др оли es ee 
рулу M | oT +h | а | k S +4 pu = - роуу!Т! = upu'v' А 
ww È (ор =) а ( ay dy ^ Рер up ) (А9) 


It is convenient, &t this point, to introduce the concepts of eddy 
viscosity, eddy thermal conductivity, and turbulent Prandtl number. The 
eddy viscosity апа eddy thermal conductivity аге defined as: 


= срам. ( A10) 


and 


ey = - РЕ (A11) 
(a 
. Ху 


The usual definition of Prandtl number is нер/Е, where the fluid 
properties are based on the molecular transport of momentum and energy. 
In this analysis there is also employed the analogous quantity Mil 
in which the properties, in effect, are based on the turbulent transport 
of momentum and energy. This quantity is called the turbulent Prandtl 
number and is given the symbol Pr ,. 


When equations (A10) and (All) are introduced into equations (A8) 
and (A9) there results 


d d d 
РУ Е = ay С + Ем) аз | (А12) 


2 

ше: а € Є | 

WW p р М 
dy 2 ау Pr Pry / dy dy 


(A13) 
Note that the bars representing time-average quantities have been dropped 


in these equations as all the terms are mean values, the fluctuating 
terms no longer appearing. 
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І.- SKIN-FRICTION RELATIONSHIP FOR Pr = Pry = 1 


Temperature-Velocity Relationship 


For the case where Pr = Pr, = 1, equation (A13) becomes 


2 
oy È (opt + Е )- E | (uve = (ет È) (Аз) 


Equation (А18) has the identical form of equation (A12). Thus, if u 
is & solution of equation (A12), then 


2 | 
ері + 25 ал + b (А15) 


is а solution of (Al!) since it is linear. The terms а end b are 
constants to be evaluated at the boundary conditions: 


Т= т at u = чу (outer edge) 
T= Ty at u = О (surface) 


When these boundary conditions are employed, equation (A15) becomes 


T ro этэр 
== 1 + Bt =- A (A16) 
м 14 


where 


ГОРЕ те» (А17) 


H 
+ 
5 


Валы UT (A18) 


Ч = — (А19) 


Because of equation (A16), it is only necessary to integrate the 
momentum equation (А12) to obtain both the velocity and temperature 
distribution in the boundary layer. To perform this integration, it is 
assumed that the boundary layer is divided into two parts: a sublayer, 
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adjacent to the surface, where momentum transfer is by molecular means 
alone; and an outer portion of the boundary layer, where momentum 
transfer is by eddying motion alone. 


Veloeity Distribution in the Sublayer 


In the sublayer, equation (A12) is written as 


11 
cr ee е А20 
Ри dy ay H 3 ) ( ) 


Since Өү is not a function of у, equation (А20) сап be integrated 
directly to yield 


рути = и B constant (А21) 


The constant in equation (A21) can be evaluated by employing the boundary 
conditions 


du 
= 0 = О = u -- А22 
Я я у-н ® (422) 
Equation (A21) becomes 
du 


Integration of equation (A23) yields 


u 

ya ЧЕ 2и (ләһ) 
о оъ + TW 

It was shown in the Discussion that the extent of the sublayer under the 

conditions of blowing is not well known. Іп View of this, extreme rigor 

in the integration of the right member of equation (A24) is not mandatory. 

For simplicity, then, let u = py and the resulting integration yields 


mo 
SF Ч : 
“> \ то ү. ер 

yt = < in [1 + tM — (А25) 
Ра 141 СР 


Ї! 


М 
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where 
5 С 
PiU, 2 y 
: ше ра (А26) 
ut = — (a27) 
СР 
JE 
Ні e (A28) 
My МТ 


Veloeity Distribution in the Outer Turbulent Portion 


In the outer turbulent portion, equation (А12) is written ав 


+ ш.а ал 
Ри dy ^ ау (s, ду (429) 


By employing Prandti's mixing-length concepts Ем = p?? Є and using 
1 = Ky, equation (A29) can be rewritten in the form 


2 


Because Qv, is not a function of у, equation (A30) can be integrated 
directly to yield 


2 
риа = py? gu | + constant (A31) 


At the interface of the sublayer and outer turbulent portion, it is 

x required that the velocity on the laminar side match the velocity on the 
turbulent side. In addition, а matching of the laminar and turbulent 
shears is also required. By comparison of corresponding terms of equa- 

x tions (A21) and (A31), it is apparent that the matching of the velocities 
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and shears requires the constants to be the same. Thus, equation (A31) 
can be rewritten as 


2 
(du 
Рима + Ty = ркеу2 au) (А32) 
Integration of equation (A32) results in 
= Уа exp [= (A33) 
E A Phy + Ty 


It will be convenient later to represent the exponential term in 
equation (A33) by the symbol g, thus 


У = Ya 6 (434) 
Because the primary concern of this analysis is the determination of 
Skin friction and heat transfer at the surface, it is not necessary to 
perform the integration indicated in equation (A33) at this stage of the 
analysis. 
Substitution of Velocity Distribution Into 
von Kármán Momentum Equation 


The von Karman momentum-integral equation for the case with surface 
blowing at constant pressure can be written as 


5 
Ty + PyVy = = | pu (u-u) dy (&35) 


Equation (A35) ean be made dimensionless by dividing through by Ра ша 9 4 
It is also convenient to replace the variable of integration у by ù. 
The resulting equation is 


f+ рМ - 4 [(g)so-w (87) ва (A36) 


The integral in the term on the right represents the momentum thickness 6. 
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To evaluate the momentum thickness rigorously in equation (A36), it 
would be necessary to divide it into two parts. The first part, with 
the lower bound аз О and the upper bound as ft, would represent the 
contribution of the sublayer. The second part, with the lower bound ав 
Ча and the upper bound as 1, would represent the contribution of the 
outer turbulent portion. In the first part (ay/at) would be evaluated 
from equation (А23), while in the second рагі it would be evaluated from 
equation (A33). For the incompressible case without blowing, these 
integrations can be performed very simply in closed form. 1% is found 
that this two-part system can be approximated very closely (less than 
l-percent difference in ce(Rg) curves) by the integral of the second part 
in which the lower bound is set equal to zero. It would not be expected 
that the introduction of compressibility and blowing should markedly 
&lter this behavior; consequently, the simpler one-part system will be 
used in the present analysis to determine 0. 


When (ау/айй) is evaluated from equation (A33), the expression for 
momentum thickness becomes 


1 3/2 
Ө = Kya "n #3 а (A37) 
Pus, St 
Оз U1 2 
noting that 
1/2 
а к( ай 
& = exp | T = (A38) 
D4U4 ча 2 


Equation (A37) cannot be solved in closed form; however, if the condition 
of compressibility is removed, a closed form solution for @ can be 
obtained easily. It is found from the solution for the incompressible 
case that the first term resulting from an integration by parts yields 
results which have the same form as the results of the complete integra- 
tion. Since we are primarily concerned with the form of the solution, 
absolute magnitudes can be adjusted by slight changes in ug and K which 
are arbitrary constants; only the first term resulting from an integration 
by parts of the compressible case will be used in this analysis. 


To perform the integration of equation (A37), it is necessary to 
replace WY as the variable of integration by в. From equation (A38) 


Аш 
Оз. 
V. С 


К 
Ча ЕЕ du | (439) 


б 
WW £ 
А pu. + 5 
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Substitution of equation (A39) into equation (A37) yields 


g(1) 


- Ya ] (&)* (14) ав (Аю) 


g(o) 


Because the pressure &cross the boundary is &ssumed to be constant, the 
perfect gas law and equation (A16) indicate that 


Р. = (т) = (а (1-80-А208)77 (ALL) 


Equation (АЮ) then becomes 


When equation (АН?) is integrated by parts, there results 


U=1 
(2 1) 1). 
= Уа NNUS А 
1-Bü-A^ 
ü-o 


From an examination of equation (A38) it is apparent that both g(1) 
and #(0) are finite; consequently, the first term on the right of 
equation (A43) vanishes at the bounds. ТР g is again substituted for 
К as the variable of integration, equation (A43) becomes 


_ Г. 1с204(А5-ВЛР au АҢ 
Ге “(1-8ВХ-Ад Су (А03) 


(n ) нө AM ЖӘЕ 
Ji g(1 ud 
T a Tu Ї [1 -eie-(A7 -B ра 0141, = 2 ds (All) 
5 ё(0) (148-28 ) 9/2 
Again, integrating by parts results in 
nG r Ч=ь 
у. AM ка 
я 3 П-о%ь(л2-ву бізі 1+5 = " 
(1481-АВТ ) 3/2 
u-o 


+ additional terms 
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When the additional terms are omitted and the limits of integration are 
introduced 


o NEA се Бү? 
" (2) ai) АБЕ и 
т ore | ae үй сэг эр (0)/-8 (АХ6) 
К (1-8-А2 575 e 
For practical cases g(1) >> (О); therefore, 
1/2 
5 Ja, Ty D + бі u и 


T ; 
But = = (1+В-А2), во that 
W 


oz% / + + (ри) g(1) (А8) 


Employing equation (A26) 


(Ako) 


(A50) 


in 


Before integrating the right member of (A50), it is preferable to obtain 
the relationship of ср in terms of the length Reynolds number. 
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Equation (A36) can be rewritten as 


when equation (A49) is employed. On integration, neglecting the effect 
of laminar sublayer | 


^ fe ит, 
i t "Eg: eig 80) ) 


с 
= > со (2. Li) 
Integr&tion by parts yields йн 
z (x) 
(Py w/ 94,1) 5 
x (@ғ/2) (1) 
Вх = + additional terms (A53) 
ди + (Puvw/ 0394) Cp 
o e 


For large values of Вх the first term predominates. Also g(1) is 
finite; therefore, 


T 
у ё(3) 
"MCN ый 
Ср ( Ср Е 
„== 12254 Lo 
2 \2 руп, 
or | "ON о «(Б ^а 
К(-- du 
e C D, V. 
“ЖУЛТ + Е 0m 
Ja 019; РИ W m: 
Piu ^ 2 


It is apparent that the right members of equations (A50) and (A55) are 
identical. This does not mean that the left members are equal, but just 
that the K and uj + шаа Ya) used in each of the equations are different 
(see Шоо: 7 
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Ihe integral appearing in equations (A50) and (A55) can be written 
as 


1/2 al 2; 
а в 
Тұ Puy Cp | 
Y. TI += J 14 BY- ABE 

1—1 


when equation (A16) is employed. The integral (A56) can be rewritten as 


т “1/2 
x (г 1 
T 
— I TERRE: NR (A57) 
a / Мих (o) со) (8-3) 
Q4 Чл 


where 
В - М B^4hA? 
a= тла (А58) 
и В+ м BA? (A 
В T од2 59) 
(сь/2) 
а = ——————— (A60) 
со 
010; 
Letting 
sin? Ф = с (в-#) ` (A61) 
then du =- = віп Ф сов Ф аф and the integral (A5T), on transforming 


the variable of integration from Ù to Ф, becomes 


| 2K (6) pe cos ф аф Tm 


pu. Фа) (вс) - 5 sin? "T (в-а) - $ sin? 3 


If c is set equal to 45 then integral (A69) becomes 
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т (ау 
Ty Ф(3) E: 7" 
--- — 1шш-ш-- 3) 
eww po (ŭa) /1 sin? Ф 
- ріс 
Ра Uy 


Since Legéndre's standard form of the elliptic integral of the first kind 
is 


Ф 
F(x, Ф) = f (лб) 
O l-k? gin? 9 


Integral (A63) can be rewritten as 


ми най 
(3 нь oso] -r е 90] } oe» 


Pu Vl HE 108 
бізі 
where 
km = 555 <1 
and | (A66) 
sin? Фа (1) = px 


If c in equation (A62) is set equal to 1/8-а, integral (A63) becomes 


2K D 
— few. то L ғы, Pn (Ma) | Ани ЕН 01] | (A67) 


A | 
22 
where 
_ B-a 
i= El 
( A68) 
sin? о, (Ф) = Bou 


NACA TN 3311 33 


II.- SKIN FRICTION - HEAT TRANSFER RELATIONSHIP FOR Pr / 1. 


This part of the analysis begins with the integration of equations 
(A19) and (A13), resulting in 


руу = (ш- Ем) = + сі (A69) 


and 


ша ©) - eB) SE orn ща 


where сі апа co аге arbitrary constants. At the surface у = 0, we 
have 


T = Ty 
du 
W = Нау 
2-2 БУР он 
aw y," Pr dy Ч 


At the surface, equations (A69) and (ATO) thus become 
O = т, + са (AT1) 
Риети = - dw + Са (А72) 


When the arbitrary constants evaluated from equations (А71) and (ATO) are 
substituted into equations (A69) and (ATO), there results 


du 
Ру Ut Ty = (Н+ ем) ау (АТЗ) 


u? T = 
гүүн | onl) +S] - а= E+ Эи 
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Dividing the terms of equation (АТ) by the terms of equation (АТЗ) 
yields 


| : цо. ем 
РУ | гу (7-2) i x | “н (ве + Pry / dep? 


^u АТ 
Puyu + Ty (H+ ey) du ыг 


It is convenient in the processes which follow to introduce the following 
symbols: 


H „ -M РъЧа 
Pr * Pr. 
2 
wW u 
= E = - — 
S Xm Cp(T-Ty) + B 


Equation (AT5) can then be written ав 


dE. PF Be PS | : 
— = Е = (1-Р) u-u ~ А 6 
In the range о<и< а Р-Рг А 


Ùa <0<1 Р = Рг} 


Equation (А76) must then be integrated in two stages with the results 
combined at the interface where Я = Ùa. The integration, with the use 
of an integrating factor, is straightforward. After some algebraic 
manipulation, there results back in the original notation 


2 Огу: 
Ма WW 
e pos 2 / piu 


па «Уты, АР сигу RPM 
pd 2 pui 2 
If we define the heat-transfer coefficient аз 
q 
Ty - Tr - гуа 


ср ср ) 
л 2 Pu 


The expression defining the recovery factor is 


Жс (сч +9) (2-Pr)(2-Pr4) - (2-Pr) Ga xs Jg) Eu 


(ыу (2 - Рг+) (2 - Pr) 


Рг-Рг+ Pry 
QV с 
s ( Pw Р st) 
Pu 2 


eol (рва: 


p 
(ка) (2 - Pr¢)(2 - Pr) 


РЕ 
И Ср t /сұуұ : 
pd “5 pup lt, 


THEE ко VOVN 


(АТ9) 
Pri 
T C 
2) -* 
2=Pry 
„+ £) Pau Se) *(Pre-Pr) 
(480) 


GE 
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Figure 1.- Direct measurements of local skin friction in low-speed flow compared with analytical 
results based on alternative choices of arbitrary parameters. 
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Figure А.- Data of reference 9 for incompressible flow with uniform transpiration 
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Figure Т.- Data of reference 9 for incompressible flow with varying transpiration rate 
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(а) М = 0; Ty = free-stream stagnation temperature. 


Figure 8.- Effect of transpiration rate on local skin-friction coefficient. 
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(4) M = 4; Ty = free-stream static temperature. 


Figure 8.- Continued. 
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Figure 9.- Effect of transpiration rate on local ratio of Stanton number to skin-friction 
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Figure 10. 


NACA - Langley Field, Va. 


